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Abstract 



Collapse and fragmentation of primordial filamentary clouds under isotropic dis- 
sociation radiation is investigated with one- dimensional hydrodynamical calculations. 
We investigate the effect of dissociation photon on the filamentary clouds with cal- 
culating non-equilibrium chemical reactions. With the external radiation assumed 
to turn on when the filamentary cloud forms, the filamentary cloud with low initial 
density (no < 10 2 cm~ 3 ) suffers photodissociation of hydrogen molecules. In such a 
case, since main coolant is lost, temperature increases adiabatically enough to sup- 
press collapse. As a result, the filamentary cloud fragments into very massive clouds 
(~ 10 5 M Q ). On the other hand, the evolution of the filamentary clouds with high 
initial density (n > 10 2 cm~ 3 ) is hardly affected by the external radiation. This is 
because the filamentary cloud with high initial density shields itself from the exter- 
nal radiation. It is found that the external radiation increases fragment mass. This 
result is consistent with previous results with one-zone models. It is also found that 
fragment mass decreases owing to the external dissociation radiation in the case with 
sufficiently large line mass. 
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1. Introduction 

It is accepted that the density perturbations collapse and cool due to hydrogen molecules 
(H2) to form so called population III (popIII) stars (Bromm et al. 1999, 2002; Abel et al. 2000, 
2002; Yoshida et al. 2008). PopIII is expected to form in halos with > 1O 6 M (Tegmark et al. 
1997). If PopIII is massive star, it is expected to affect neighbor clouds via radiative feedbacks. 
Radiative feedbacks cause ionization and dissociation. Although ionization photon tends to be 
prevented from spreading out of halos because of large opacity of hydrogen atoms, dissociation 
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photon tends to spread out of halos (Kitayama et al. 2004). Thus, some regions are expected 
not to be ionized but to be photodissociated. We consider the filamentary clouds in such a 
region. 

Filamentary clouds are possible origin of stars. In general, non-spherical gas cloud 
tends to become sheet-like cloud, and sheet-like cloud tends to fragment into the filamentary 
clouds (Miyama et al. 1987). In numerical cosmological simulations of first star formation, 
filamentary structure is frequently seen (e.g., Abel et al. 1998; Bromm et al. 1999; Greif 
et al. 2006). Recently, many filamentary structures have been found through the Herschel 
Gould Belt Survey (Andre et al. 2010). A filamentary cloud is possible to fragment into many 
quasi-spherical clouds (Nagasawa 1987; Inutsuka & Miyama 1997). These spherical clouds 
are expected to become stars or other astronomical objects. We investigate how this process 
exceeds when first stars form. 

There are previous works about fragmentation of the filamentary cloud (Uehara et al. 
1996 ; Nakamura & Umemura 1999, 2001, 2002; Flower 2002; Omukai & Yoshii 2003). Among 
these, Nakamura & Umemura (1999, 2001, 2002) used one- dimensional hydrodynamical calcu- 
lations and two-dimensional hydrodynamical calculations. The authors considered many cases 
with various initial density, temperature, line mass, and initial fraction of H 2 . It is found 
that fragment mass is 1 — 5OOM and has a bimodal distribution when initial H2 fraction is 
10 -3 . However, the effect of the external radiation was not considered. Hence, these works 
is applicable only to first star formation. Among works mentioned above, Omukai & Yoshii 
(2003) considered the external dissociation radiation. Their work is applicable to formation of 
second-generation stars. They calculated the thermal evolution of the filamentary cloud under 
the isotropic external radiation with one-zone model assuming free-fall. They assumed that 
the filamentary cloud fragments when its density becomes 100 times higher than the loitering 
point. Under this assumption, Omukai & Yoshii (2003) concluded that the effect of the external 
dissociation radiation decreases fragment mass. 

In Bessho & Tsuribe (2012) (hereafter Paper I), we investigated collapse and fragmenta- 
tion of a filamentary clouds under the isotropic external radiation using one-zone models. We 
assumed that the external radiation turns on when the filamentary cloud forms. By taking into 
accounts pressure effect explicitly, we found that the filamentary clouds with low initial density 
(no < 10 2 cm -3 ) suffer photodissociation and fragment into very massive clouds (~ 1O 4_5 M ). 
It was found that the effect of the external radiation increases fragment mass. The evolution 
of the filamentary clouds with high initial density (no > 10 2 cm~ 3 ) or with sufficiently large 
line mass is not affected by the external radiation owing to self-shielding. In Paper I, at first, 
we assumed the uniform filamentary clouds with homologous collapse. However, in realistic 
situations, filamentary clouds are expected to collapse in run-away fashion. Hence, in Paper I, 
we also introduced "rarefied filament model" as an improved model which includes the effect of 
run-away collapse partly. In rarefied filament model, line mass of collapsing core decreases as 
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rarefaction wave propagates from the cloud surface. As a result of the rarefied filament model, 
we found that the effect of the external radiation increases fragment mass. However, this result 
is apparently inconsistent with Omukai & Yoshii (2003). 

The purpose of the present paper is to investigate whether or not fragment mass in- 
creases owing to the effect of the external radiation and to clarify the reason why our result 
is inconsistent with Omukai & Yoshii (2003). For this purpose, we extend our previous in- 
vestigation using the one-dimensional model which includes the full characteristics of run-away 
collapse and have more realistic results. Furthermore, in this paper, we also consider the further 
evolution of each fragment. 

We assume that the external radiation is isotropic. Intensity of the external dissociation 
radiation is set according to distribution of intensity of the dissociation radiation at redshift 
z ~ 10 (Dijkstra et al. 2008). The intensity which we consider is moderate (§2.2), and we 
consider the situation where the dissociation radiation originates from halos out of a halo 
including the filamentary cloud. 

In Paper I, we concentrated on the case where the external radiation turns on when 
the filamentary cloud forms in > 10cm~ 3 ). To clarify the reason of the apparent difference of 
the conclusions between us and Omukai & Yoshii (2003), we also consider the case where the 
external radiation turns on at lower density (n = 0.1cm -3 ) as in Omukai & Yoshii (2003). 

In §2, we describe our model for the filamentary clouds. We present numerical results in 
§3. In §4, we investigate the reason for the difference between us and Omukai & Yoshii (2003). 
§5 is devoted to conclusions and discussion. 

2. Model for the filamentary clouds 

2.1. Basic equations 

We assume the axisymmetric filamentary cloud. We do not consider dark matter for 
simplicity. This simplicity gives us a good approximation in the case with high initial den- 
sity. In the case with low initial density, the effect of dark matter gravity is underestimated. 
Hydrodynamical equation of motion in Lagrangian form is given by 

2L = -™ (!) 

Dt r or K J 

where v is velocity in the cylindrical radial direction, G is gravitational constant, I is line mass 
within cylindrical radius r, and P is pressure for ideal gas given by 

P = nk B T, (2) 

with number density n, temperature T, and Boltsmann constant ks- 
We also solve energy equation given by 

dt~ dtp p U 



where p is density and u is thermal energy per unit mass, 
1 k B T 

u = 4 

7 ad - 1 /im H 

with adiabatic index 7 ac j, mean molecular weight fi, and mass of a hydrogen atom ran. The 
symbol A net in equation (3) is cooling rate per unit volume including lines of H, lines of H 2 , lines 
of HD, and chemical heating/cooling. Since continuum processes hardly change the evolution 
of the filamentary clouds (Paper I), we neglect them. As for lines, we estimate cooling rate 
from detailed balance of population of energy levels. The escape probability for emission by 
the transition between level i and j is given by 
1 — ft - ry 

h = — —> (5) 

with assuming the velocity profile v r (r) oc r (Castor 1970). The cooling rate is multiplied by 
this escape probability. Optical depth, r^ , is given by 

rRout 

= / Kij(r') dr' 



Rout flVjA 

— (rijBii — riiBi~\dr' (6) 

where R out is radius of the outer boundary, «y is opacity for lines, hv^ is energy difference 
between levels % and j, rij (rij) is the level population at level i (j), Bij and Bji are the Einstein 



-B-coefficients, and Au^ = j ' c^jlksT '/ 'A* m H is the thermal Doppler width of transition line 
i-tj. 

We consider non-equilibrium chemical reactions by solving following equation for each 
fluid element : 

777 = £ kijkfjfkn + k *jfh ( 7 ) 
ai j, k j 

where kijk and kij are reaction rates for formation and destruction of speices i, and fi is fraction 
of species i. We consider fourteen species : H, H + , H~, H 2 , H^, He, He + , He ++ , D, D + , D _ , 
HD, HD + , and e™. We consider 26 chemical reactions concerned with H and He taken from 
Nakamura & Umemura (2001), photodissociation of H 2 (equation 10), and 18 chemical reactions 
concerned with H and D taken from Nakamura & Umemura (2002). We solve equation (7) with 
implicit integrator. 

We solve equations (l)-(7) using 200 spatial meshes in cylindrical radial direction. As 
for initial interval of meshes, we set Arj + i = l.OlAr,;, where Arj = rj+i — 7*. Hence, spatial 
resolution is much better in the central region. The mesh size Ar^ is checked to be shorter than 
1/4 of local Jeans length 1 at all times (Truelove et al. 1997). We set outer boundary to be 



The exact definition of Jeans length is given by 
\GjjLmnn 



/ \ 1/2 

Aj= 7^ c s- (8) 
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Fig. 1. Error of density distribution of isothermal equilibrium filamentary 
cloud (T = 300K) for the numerical results with 100 meshes and 200 
meshes. The symbol p ana is analytic solution, and p com is numerical results. 



10i?o as a sufficiently large value, where Rq is the effective radius given by 

Rq ■ 



2fk B T 



(9) 



7rG/i 2 m^n ' 

with line mass parameter / (equation 17), initial temperature To, and initial number density 
at the center uq. As for the outer boundary condition, we assume that the external pressure is 
zero. 

We solve equation (1) in the second-order-accurate finite-difference scheme with the 
artificial viscosity (Richtmyer & Morton 1967). We calculate in the same way as Thoul & 
Weinberg (1995), except for cylindrical geometry. To check accuracy of our code, first, we 
calculate the density distribution of isothermal equilibrium filamentary cloud (T = 300K) with 
100 and 200 meshes, and fixed time step. Drag term —2v(r)/dt is added to equation of motion 
and is eliminated after the step number reaches 125 for case with 100 meshes and 250 for case 
with 200 meshes. After 250 steps (100 meshes) and 500 steps (200 meshes), we have error of 
density distribution shown in figure 1. Error for case with 100 meshes is 4 times of error for 
case with 200 meshes. Hence, our code is second-order-accurate in space. 

Second, we check temporal accuracy of our code in time. We calculate collapse of 
pressure-less uniform filamentary cloud in free-fall state. Analytic solution is given by 



v / - lo e F W _ 2 

e x dx 



o 



nGp(0)t, 



(10) 



where F(t) = r(i)/r(0), p(0) is the initial density. Time step is set to be At = 10~H s (t = 0) 
and 2 At. Calculation is continued until the density becomes 100 times of the initial density. 
In the case with At, F at the end of calculation is 0.0913290. In the case with 2 At, F at the 
end of calculation is 0.09132307. Analytic solution predicts F = 0.9132284. Error with 2At is 
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4 times larger than At. Hence, our code has second-order-accuracy in time. 

2.2. External radiation 

We assume that the external dissociation radiation is isotropic. Intensity of the external 
radiation is set according to Dijkstra et al. (2008). Dijkstra et al. (2008) calculated the 
probability distribution of mean intensity of the dissociation radiation at redshift z ~ 10 by 
estimating mean intensity of the dissociation photon from the surrounding halos to a single 
halo. In this paper, we use mean intensities whose probability is 0.4 and 0.06 as in Paper I. 
The external radiation is assumed to be thermal radiation of 12OM star, and we determine 
surface temperature (T sur = 95719K) according to Schaerer (2002). We assume that the external 
radiation turns on when the filamentary clouds forms in §3. In §4, we change the density when 
the external radiation turns on. 

We calculate the photodissociation reaction of H 2 , 

H 2 + 7 ^H 2 *^2H, (11) 

(solomon process) where 7 is photon with 12.4eV and H 2 * is excited state of H 2 . The reaction 
rate is given by 

fc 2step = 1-4 X 10 9 /shJext S"\ (12) 

where J ext is mean intensity of the external radiation at the surface of the filamentary cloud 
and / s h is self- shielding function, 

^ = min [ 1 ' (10^) 'I' (13) 

where Nr 2 is column density of H 2 (Draine & Bertoldi 1996) estimated as 

N U2 (r) = J r R ° Ut nn 2 (r')alr'- (14) 

2. 3. Fragmentation of the filamentary cloud 

There are two important timescales during the collapse of filamentary clouds. One 
is timescale for density evolution, t^ yn = p c /p c , and the other is timescale for fragmentation, 
i frag = 5.17/ \J2ixpcG (Nagasawa 1987). The latter is the timescale during which the fastest 
growing mode of perturbation grows to non-linear. If the fastest growing mode has time enough 
to grow to non-linear before the fastest growing mode changes, the filamentary cloud is expected 
to fragment. Thus, we assume that the filamentary clouds fragment when tdyn > ifrag- Once 
condition for fragmentation is satisfied, it has been satisfied after that. 

We estimate fragment mass by integrating region with the density higher than 10% of 
central density n c , 

rr(n=0.1n c ) 

M frag = A frag / 2nr'pdr', (15) 

J 

where 



g 0.288v/4vrG>m H n c v ; 

is the wave length of the fastest growing mode for equilibrium filamentary clouds (Nagasawa 
1987) and c s is sound speed. Since Af rag oc p -1 ^ 2 , fragment mass is smaller when the filamentary 
cloud reaches higher density before fragmentation. Since interval of integration in equation 
(14) is approximately Jeans length, Jeans mass at fragmentation is close to fragment mass (see 
§3.3). 

2.4- Parameters and initial conditions 

In this paper, we treat three physical quantities as parameters. First is initial density 
no, second is normalized mean intensity of the external radiation, 

t Jfiv=13.6eV,ext »\ 

and third is line mass parameter 2 , 

/= 2k B T R °- (18) 
Line mass parameter is important in the view point of dynamical evolution. Initial density no 
is important for thermal evolution. With respect to photodissociation, mean intensity J21 and 
initial density n$ are important. 

We consider cases with log 10 ra = 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, and 6 for n and 
/ = 1.5, 2, 2.5, 3, 3.5, 4, 5.5, and 6 for /. For J 2 i, we consider J 2 i = 0, 1, 6.5, and 10. According 
to Dijkstra et al. (2008), the case with J21 = 1 represents the weak external radiation case, 
J21 = 6.5 is mean intensity with the highest probability (0.4), and J21 = 10 represents the strong 
radiation case whose probability is 0.06. 

We assume initial density distribution as 

n(r)=n (l + ^J (19) 

with Rq given in equation (8). With / = 1, equation (18) represents equilibrium density distri- 
bution for the isothermal filamentary cloud (Ostriker 1964). In this paper, we concentrate on 
collapsing filamentary cloud with / > 1. Initial velocity distribution is assumed to be 

«(r) = r=r=J- ( 2 °) 

Ro + y/R 2 + r 2 

Equation (19) indicates that infall velocity is proportion to radius at r — > and is constant 
(~ —c s ) at r — > 00. Although the actual initial velocity may be different from equation (19) 
depending on the detail of dynamical evolution of the filamentary cloud formation, results in 
this paper will not change qualitatively unless initial velocity is much faster than a few times 
of equation (19). 

2 If a filamentary cloud forms as a result of fragmentation of the fastest growing mode in a sheet-like cloud, 
the typical value of / is 2 (Miyama et al. 1987). 
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Fig. 2. Evolution of the temperature (a), /h 2 and /hd X 10 4 (b), the heating and cooling rate 
(c), and idyn and if rag (d), respectively, as a function of the central density for the case with 
(/, no, J 21) = (1.5, 10cm -3 , 0). In the diagram (c), "adiabatic" denotes the adiabatic heating, "H 2 " does 
the H2 line cooling, "HD" does the HD line cooling, and "chemical" does the chemical heating or cooling. 

As initial temperature, we adopt T = 300K, assuming that the filamentary cloud forms 
after the cloud undergoes H 2 cooling. We also adopt /h 2 = 10~ 4 and / e = 10~ 4 . As for the 
value of /h 2 , we refer the result in Paper I. Fraction of electron is set in order not to change 
/h 2 artificially via H~ channel. Initial fraction of proton is set to be f p = 10 -4 for charge 
conservation. We assume [H]/[D] = 4 x 10~ 5 , which is consistent with observations of the 
deuterium Lya feature (e.g., O'Meara et al. 2001). Initial fraction of the others is set to be 
zero. 

3. Results of one-dimensional hydrodynamical calculations 

3.1. Cases without the external radiation 

To investigate the effect of the external radiation, at first, we show the results of the 
cases without the external radiation. 

3.1.1. Low density filamentary cloud with small line mass 

First, we show the results for the case with low initial density and small line mass, 
(/, no, J21) = (1.5, 10cm -3 , 0) (figure 2). In the early stage of collapse, H 2 cooling dominates 
adiabatic heating a little, and temperature decreases. After the density reaches n ~ 10 3 cm -3 , 
adiabatic heating dominates H2 cooling. In n > n cr it ~ 10 4 cm~ 3 , H2 cooling rate is proportional 
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n (cm 3 ) 

Fig. 3. Same as figure 2, but (/,n , J 21) = (6, 10cm~ 3 ,0). 

to n, while it is proportional to n 2 in n < n crit . Hence, since cooling time becomes constant 
in n > n cr it, dynamical time becomes constant and longer than fragmentation time. When 
n reaches ~2x 10 5 cm -3 , condition for fragmentation is satisfied with Mf rag ~ 122OM . To 
ensure that fragmentation occurs, we continue to calculate the evolution until free-fall time has 
past after the condition for fragmentation is first satisfied. Once condition for fragmentation is 
satisfied, it has been satisfied. 

3.1.2. Low density filamentary cloud with large line mass 

Next, we show the result for the case with low initial density and large line mass, 
(f,no, J21) = (6, 10cm -3 , 0) (figure 3). Owing to larger line mass, the filamentary cloud collapses 
to higher density than figure 2. Collapse continues up to high density (n ~ 10 6 cm -3 ), and 
the filamentary cloud fragments. Fragment mass is ~ 49OM . Until fragmentation, adiabatic 
heating and H2 cooling balance, and temperature is approximately constant (T ~ 350K). 

3.1.3. High density filamentary cloud with small line mass 

As the final example, we show the result for the case with high initial density and small 
line mass, (f,n ,J 2 i) = (1.5, 10 6 cm~ 3 , 0) (figure 4). In the early stage of collapse, Adiabatic 
heating dominates H 2 cooling and temperature increases. During collapse, H 2 cooling rate 
increases and approximately balances with adiabatic heating rate at n ~ 3 x 10 6 cm~ 3 . After 
temperature decreases a little, fragmentation condition is satisfied at n~ 10 7 cm -3 since collapse 
is suppressed owing to high temperature. Fragment mass is ~ 37OM . 

In summary, in the cases without the external radiation, the filamentary cloud undergoes 
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Fig. 4. Same as figure 2, but (f,n , J 21) = (1.5,10 6 cm 3 ,0). 



approximately isothermal states and fragments. This feature comes from the fact that H 2 
cooling and adiabatic heating compete each other. The results in this subsection are similar 
to previous works (Nakamura & Umemura 2001, 2002). For the parameter set same as figure 
2, in Paper I, fragment mass was 23M Q in the uniform model and 3500M Q in rarefied filament 
model. Fragment mass 122OM in the one-dimensional model is close to the results of the 
rarefied filament model. This result indicates that the effect of run-away collapse is important 
to estimate fragment mass. 

3.2. Cases with the external radiation 

In this subsection, using the same parameters as figure 2, 3, and 4, we investigate how 
the external radiation changes the thermal evolution and fragment mass of the filamentary 
cloud. 

3.2.1. Low density filamentary cloud with small line mass and strong radiation 

First, we show the result for the case with low initial density, small line mass, and 
strong external radiation, (f,no, J21) — (1.5, 10cm~ 3 , 10) (figure 5), where the external radiation 
is added to the case of figure 2. This case would be affected by the external radiation because 
of low density. Most of H2 are photodissociated in the early stage of collapse, and temperature 
increases adiabatically. The filamentary cloud fragments into very massive fragments (~ 2.4 x 
10 5 M Q ) at n ~ 30cm~ 3 . This result demonstrates that in the case with low initial density 
and small line mass the external radiation increases fragment mass and changes the thermal 
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Fig. 5. Same as figure 2, but (f,n , J 2 i) = (1.5, 10cm -3 , 10). 

evolution. 

3.2.2. Low density filamentary cloud with large line mass and strong radiation 

Next, we show the result for the case with low initial density, large line mass, and strong 
external radiation, (f,n , J 2 i) = (6, 10cm -3 , 10) (figure 6), where the external radiation is added 
to the case of figure 3. This case is also excepted to be affected by the external radiation because 
of low density. However, the filamentary cloud may collapse up to higher density than figure 
5 because of large line mass. In figure 6, it is seen that most of H2 is photodissociated in the 
early stage of collapse and temperature increases adiabatically as in figure 5. However, since 
the filamentary cloud is more massive than in figure 5, stronger gravity and large inertia help 
collapse. Fragmentation does not occur during early adiabatic phase, and collapse continues 
until the density becomes higher than in figure 5. At n~ 10 2 cm -3 , H2 starts to form and shields 
itself from the external radiation. The filamentary cloud starts to cool owing to H2 cooling. 
After n ~ 10 3 cm -3 , Since H 2 cooling balances with adiabatic heating, temperature becomes 
nearly constant (T ~ 400K). The filamentary cloud fragments into clouds with ~ 590M Q at 
n ~ 10 6 cm~ 3 . Fragment mass and density at fragmentation are similar to the case without the 
external radiation (figure 3). In the case with large line mass (/ = 6), it is found that fragment 
mass is hardly affected by the external radiation although the evolution of temperature is 
affected by the external radiation in the early stage of collapse. 
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Fig. 6. Same as figure 2, but (f,n , J 21) = (6,10cm 3 ,10). 
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Fig. 7. Same as figure 2, but (f,n ,J 2 i) = (1.5, 10 6 cm- 3 , 10). 
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3.2.3. High density filamentary cloud with small line mass and strong radiation 

Finally, we show the result for the case with high initial density, small line mass, and 
strong external radiation, (f,no, J21) = (1.5,10 6 cm~ 3 ,10) (figure 7), where the external radiation 
is added to the case of figure 4. This case may not be affected by the external radiation 
because of high density. Since initial density is high enough to shield the filamentary cloud 
from the external radiation, H2 near the center of cloud is not photodissociated. The evolution 
of temperature at the center is hardly affected by the external radiation and is similar to figure 
4. The filamentary cloud fragments into clouds with ~ 400M© at n ~ 10 7 cm~ 3 . It is found 
that the effect of the external radiation is not important in the case with high initial density 
(n = 10 6 cm~ 3 ). 

In summary, in the case with low initial density (uq < 10 2 cm~ 3 ), the filamentary cloud 
suffers photodissociation in the early stage of collapse. In such case, temperature increases 
adiabatically. The filamentary cloud with small line mass (/ = 1.5) fragments during adiabatic 
phase. On the other hand, the filamentary cloud with large line mass (/ = 6) does not fragment 
during adiabatic phase and collapses with shielding itself from the external radiation. In this 
case, fragment mass is hardly affected by the external radiation. In the case with high initial 
density (no = 10 6 cm~ 3 ), the thermal evolution of the filamentary cloud is hardly affected by 
the external dissociation radiation. 

One-zone model predicts fragment mass different from 2.4 x 10 5 M Q which one- 
dimensional model predicts. For example, for the parameter set same as figure 5, in Paper I, 
the uniform model predicted 1.5 x 1O 5 M and rarefied model predicted 2.7 x 10 4 M Q . Difference 
between one-zone models and one-dimensional model originates from difference of dynamical 
equation (virial equation in one-zone model and hydrodynamical equation of motion in one- 
dimensional model). In one-dimensional model, collapse is run-away collapse, and fragmenta- 
tion condition is satisfied at lower density since free-fall time balances with sound crossing time 
in the central dense region. Hence, in one-dimensional model, fragment mass becomes larger 
than in one-zone model. 

3. 3. Property of the filamentary cloud at fragmentation 

In this subsection, we show the profile of physical quantities (density, temperature, in- 
fall velocity, and ratio of pressure gradient to gravitational force) at fragmentation. We focus 
on density profile and investigate whether or not the universal profile at fragmentation exists. 
Moreover, we compare fragment mass (equation 14) with Jeans mass estimated with the central 
density and temperature. 
3.3.1. Case without the external radiation 

We show the profiles of density, temperature, infall velocity, and ratio of pressure gradient 
to gravitational force at fragmentation in the case with (/, n , J21) = (1.5, 10cm -3 , 0) (figure 
8). In diagram (a) in figure 8, it is seen that dense central region within Jeans length (Aj ~ 
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r (cm) r (cm) 

Fig. 8. Profile of density (a), temperature (b), in-fall velocity profile (c), and ratio of pressure gra- 
dient to gravity (d) at fragmentation of the filamentary cloud with (/, no, J21) = (1.5, 10cm -3 , 0). 

2.3 x 10 18 cm) has uniform density, and the density profile in the outer envelope is proportional to 
r -4 . This density profile is similar to that of equilibrium solution for the isothermal filamentary 
cloud (Ostriker 1967). However, between r = 10 19 cm and r = 10 20 cm, slope of the density profile 
is sallower than r~ 4 . Temperature is highest outside r coo i ~2x 10 18 cm where t coo i = tg, and 
pressure gradient force is stronger than gravity force outside r coo i. Hence, matter is pushed 
outward. Velocity profile is in proportion to radius in the central dense region and is constant 
larger than sound speed in the outer envelope. Ratio of pressure gradient to gravity is nearly 1 
(~ 1.01) inside r coo i. In diagram (b), drop of temperature at the surface is seen. Since we assume 
that the external pressure is zero, adiabatic cooling occurs at several meshes of the surface. 
Moreover, these meshes are pushed by inner meshes with higher pressure and fall more slowly 
than inner meshes. However, these effects do not affect the central region. Fragment mass is 
1220M Q which is close to Jeans mass (114OM ) estimated with central density and temperature. 
Since tdyn is about 6 times of ts at the center when the filamentary cloud fragments, pressure 
gradient force is important to calculate the further evolution of fragments. Further evolution 
of each fragment is shown in §4. 
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Fig. 9. Same as figure 8, but (f,no, J21) = (1.5,10cm 3 ,10). 



3.3.2. Case with the external radiation 

We show profiles of the same quantities as figure 8 for the case with (/, no, J 21) = 
(1.5, 10cm -3 , 10) in figure 9. In figure 9, it is seen that except for temperature, profiles of 
physical quantities are similar to figure 8. Most of H 2 is photodissociated, and the filamentary 
cloud loses the ability to cool. Hence, temperature is higher in the central dense region than 
in the outer envelope. Ratio of pressure gradient to gravity is larger than 1 (~ 1.2). Fragment 
mass is 2.4 x 1O 5 M which is close to Jeans mass (3.0 x 1O 5 M ) estimated with central density 
and temperature. In the case with the external radiation, since tdyn is about 5 times of tg at the 
center when the filamentary cloud fragments, pressure gradient force is important to calculate 
the further evolution of fragments (see §4). In figure 10, the density profiles in figure 8 and 9 
are simultaneously plotted, and each profile is found to be similar to each other. The profiles 
at r/Xj < 0.5 are similar to the profile of isothermal filamentary cloud in equilibrium state. 

3-4- Fragment mass 

We show how much the external dissociation radiation changes fragment mass. Figure 
1 1 shows the fragment mass for all the parameters in n — / plane using contours. Results for 
the cases with J21 = 0, 1, 6.5, and 10 are presented in different diagrams. In the case with 
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figure 8 




Fig. 10. Normalized density profile for figure 8 (solid line) and figure 9 (dashed 
line). The initial density profile in the case with / = 1 is also shown (dot- 

ted line). The symbol n c is central density and Aj is Jeans length at the center. 
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Fig. 11. The contours map for the fragment mass for the case with (a) J21 = 0, (b) 
J21 = f, (c) J21 = 6.5, and (d) J 2 i = fO. Solid lines in each diagrams represent con- 
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the external radiation, it is seen that the filamentary clouds fragment into very massive clouds 
(> 1O 5 M ) in the cases with low initial density (uq < 10 2 cm~ 3 ). Since very massive fragments 
are not seen in the case without the external radiation, this can be regarded as a result of the 
effect of the dissociation photon. This feature is similar to the result of the one- zone models in 
Paper I. Thus, formation of very massive fragments under the external radiation with moderate 
intensity can be regarded as the robust result, provided that the external radiation turns on 
when the filamentary cloud forms. 

The diagram (a) of figure 11 is similar to figure 6 of Nakamura & Umemura (2001). In 
the range of 2 — 100M o , the contours are dense. This is because the filamentary cloud becomes 
isothermal once H2 cooling becomes effective owing to three body reaction, and continues to 
collapse to high density (~ 10 13 cm~ 3 ). In such a case, fragment mass is small (~ 2 — 10M Q ). 
This feature is seen in Nakamura & Umemura (2001). 

Nakamura & Umemura (2002) concluded that there are some parameter sets where HD is 
main coolant. However, in our results, HD is found not to be important. Deuterated hydrogen 
molecules HD mainly forms from H2 3 . Since initial H2 fraction in this paper is assumed to be 
small (10 -4 ), even in the case without the external radiation, sufficient amount of HD to cool 
does not form. This result is consistent with Nakamura & Umemura (2002). In the case with 
the external radiation, since H2 is photodissociated, HD is less important. 

We show the effect of the external radiation on fragment mass quantitatively. Figure 12 
shows the similar contours but about the ratio of the fragment mass between the cases with 
and without the external radiation. In addition to figure 11, figure 12 clearly shows that the 
filamentary clouds with low initial density (n < 10 2 cm -3 ) and moderate / (< 4.5) fragment 
into more massive fragments than the case without the external radiation. This feature agrees 
with the results of the rarefied filament model in Paper I. 

It is seen that fragment mass for the case with high initial density (no > 10 2 cm -3 ) does 
not increase owing to the external dissociation radiation. This is because the initial density is 
high enough for the filamentary cloud to shield itself from the external dissociation radiation. 
In this case, the evolution is similar to the case without the external radiation. 

3.5. Self-shielding function of Wolcott- Green et al. (2011) 

Wolcott-Green et al. (2011) have recently suggested self- shielding function indicated by 
three-dimensional radiative transfer. In this subsection, we investigate dependence of fragment 
mass on the self- shielding function quantitatively. The new self-shielding function is given by, 

/ sh ,wc = + (T^p ex P[- 8 - 5 x + *)°- 5 ]> ( 22 ) 

3 Main chemical reaction of formation of HD is given by 

H 2 +D+^HD + H+. (21) 
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Fig. 12. The ratio of the fragment mass with the external radiation ((a) 
J21 = 1, (b) J21 = 6.5, and (c) J21 = 10) to that without the external radia- 
tion {J 21 — 0). The number near each solid line is the ratio of fragment mass. 



where 



and 



x = 



H 2 



5 x 10 14 cm- 2! 



1 



l2k B T 



(23) 



(24) 



10 5 cm/s 10 5 cm/sV fJ-rriH 
Since / s h,WG is 1 — 10 times larger than original / s h (equation 12). Hence, Wolcott-Green et al. 
(2011) suggested that the effect of photodissociation is actually stronger than the result with 
/ s h- Thus, the results which we have ever considered are expected to be modified quantitatively. 

To see the difference between shielding function, we calculate fragment mass using / s h,wc 
in the case with J21 = 10. Figure 13 shows the ratio of fragment mass between the cases with 
/sh,cw and with / S h. When we use / S h,GW, fragment mass increases comparing with the case 
with / s h especially for the cases with no = 10 — 10 3 cm -3 and / > 2. This feature is consistent 
with the relation, / sh , G w > /sh- 
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Fig. 13. The ratio of the fragment mass with jsh.GW to that with / s h in the case 
with J21 = fO. The number near each solid line is the ratio of fragment mass. 

4. Criterion for increase of fragment mass 

Omukai & Yoshii (2003) calculated the evolution of the filamentary cloud under the 
external dissociation radiation assuming free-fall. The authors assumed that fragmentation 
occurs at density 100 times higher than the loitering point and concluded that the effect of the 
external dissociation radiation decreases fragment mass. This conclusion apparently disagrees 
with our results in §3. In this section, we consider whether or not the external radiation 
increases fragment mass when the filamentary cloud reaches the loitering point. Initial density 
of the filamentary cloud is assumed to be very low (no = 0.1cm -3 ). As the timing when the 
external radiation turns on, we consider various cases. The investigation in this section provides 
systematic study which includes the situation of Omukai & Yoshii (2003). 

4-1. Whether the filament reaches the loitering point 

Suppose a filamentary cloud with no = 0.1cm -3 and the external radiation turns on when 
density reaches nuv- We consider the cases with nuv = 0.1cm -3 , lem -3 , and 10cm -3 . We also 
assume T = 300K and / H2 = at n . As for the evolution of the filamentary cloud, we solve 
one- dimensional hydrodynamics as in §3. 

We calculate the evolution of the filamentary cloud with various values of /. In order 
for the filamentary cloud not to fragment during adiabatic phase in the case with J 2 \ = 10, 
it is found that / is required to be larger than 30, 25, 10, and 5 for various values of nuv 
; nuv — 0.1cm -3 , 1cm -3 , 10cm -3 , and 00, respectively. Hence, we investigate how massive 
fragments are in the case with / = 30, 25, 10, and 5. Does fragment mass increase when the 
filamentary cloud reaches the loitering point? 

First, we show the case with / = 30 where the filamentary cloud reaches the loitering 
point for any nuv- Figure 14 shows the thermal evolution of the filamentary cloud with / = 30, 
n = 0.1cm -3 , J 2 i = 10, and various nuv- Thin lines indicate thermal evolution of each fragment 
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Fig. 14. Thermal evolution of the filamentary cloud with / = 30, no = 0.1cm -3 , J^x = 10, and var- 
ious iiuv where the external radiation turns on ; nuv = 0.1cm -3 (thick solid line) titjv = 1cm -3 
(thick long-dashed line), nuv = 10cm -3 (thick short-dashed line), nuv = oo (thick dash-dotted line). 
Thin lines shows the thermal evolution of fragments. Dotted lines indicate the constant Jeans masses. 




Fig. 15. Same as figure 14, but / = 5. 



(see §4.2). Fragment mass is largest (550M Q ) in the case with nuv = oo. Fragment mass is 
smaller in the case with lower nuv- As a result, the external radiation decreases fragment 
mass. All filamentary clouds fragments at the density which is roughly 1000 times higher than 
loitering point, which is qualitatively consistent with Omukai & Yoshii (2003). 

Next, we investigate whether or not the external radiation decreases fragment mass when 
the filamentary cloud does not reach the loitering point. Figure 15 shows the same figure as 
figure 14 except for the value of / : / = 5. It is found that the filamentary cloud which does not 
reach the loitering point fragments into more massive clouds than the case with nuv = oo. This 
feature agrees with our result of §3. Also in the cases with f = 25 and 10, the filamentary cloud 
fragments into more massive clouds than the case without the external radiation (nuv = oo) 
when it does not reaches loitering point (nuv < 0.1cm -3 for f = 25 and nuv < 1cm -3 for / = 10). 
Whether or not the external radiation increases fragment mass is determined by whether or 
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not line mass is small enough for the filamentary cloud not to reach the loitering point. 

It is explained as follows that fragment mass decreases under the external radiation : 
after the external radiation turns on, the filamentary cloud suffers photodissociation. The outer 
region loses more amount of H2 than the dense central region. Since temperature of the outer 
region becomes higher than the dense central region, pressure gradient becomes sallower at the 
region between the dense central region and the outer region. Then, such a region sticks to the 
dense central region owing to gravity, and line mass of the central dense region increases. As 
a result, the filamentary cloud collapses up to higher density (n ~ 10 6 ~ 7 cm~ 3 ) owing to larger 
line mass and fragments into less massive clouds than the case without the external radiation. 

4-2. Sub-fragmentation 

It is found that the external radiation increases fragment mass when the filamen- 
tary cloud fragments before the loitering point. However, each fragment may occur sub- 
fragmentation and mass of final outcome may be as small as fragment mass in the case without 
the external radiation. To clarify mass of final outcome, we investigate the further evolution of 
each fragment. 

We use one- zone model for each fragment. We assume that each fragment is spherical 
and has Mf rag (see equation 14), chemical composition same as at the center of the filamentary 
cloud, radius ro.i where n(ro.i) = 0.1n c . Infall velocity at fragmentation is proportional to 
radius, and we use infall velocity at r .i as infall velocity of each fragment. The density of 
fragment is estimated from the mass and radius as 

Pfrag = (25) 
— r 0.1 

This density, Pf rag , is different from the density of the filamentary cloud before fragmentation 
and corresponds to the one at the end of fragmentation. Temperature and fraction of each 
fragment is approximated by the value when the filamentary cloud fragments. According to 
the result in §3.3, in the view point of dynamical evolution, we should consider the effect of 
pressure gradient force. Hence, we treat pressure effect explicitly using the virial equation for 
uniform sphere. The virial equation is given by 

dv = 10 ]__GM_ 

alt 3 iimn R R 2 

where R is radius and M is mass of the cloud (see Appendix 1). Radiative transfer is treated 
by the same method as in Omukai (2001) except for we use / s h,cw instead of / s h- The same 
routine as §3 is used for chemical reactions. 

In figure 15, it is found that if each fragment occurs sub-fragmentation at the loitering 
point (n ~ 10 3 cm -3 ), mass of final outcome is ~ 1O 4 M . This is larger than fragment mass in 
the case without the external radiation. Hence, even if sub-fragmentation occurs, mass of final 
outcome increases owing to the external radiation when the filamentary cloud fragments before 
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the loitering point. This tendency is found in the other cases with / = 25 and 10. 

When the filamentary cloud reaches the loitering point, temperature of each fragment 
increases adiabatically after fragmentation {n ~ 10 6_7 cm -3 ). This is because gravity dominates 
pressure gradient after fragmentation and dynamical time becomes shorter than cooling time. 

5. Conclusions and discussions 

In this paper, collapse and fragmentation of primordial filamentary cloud was investi- 
gated using one-dimensional hydrodynamical calculations with the effect of the external disso- 
ciation radiation. Especially, the effect of run-away collapse to fragment mass is considered by 
comparing with previous results with one-zone models. Results are summarized as follows : 

• Comparing with the uniform model in Paper I, one-dimensional filament model predicts 
lower fragmentation density and larger fragment mass. This is because fragmentation 
occurs only in the central region with low virial temperature in one-dimensional model. 

• Comparing with the rarefied filament model in Paper I, one-dimensional filament model 
predicts similar fragment mass. This explains that the discrepancy between the uniform 
filament model and one-dimensional filament model mainly comes from the run-away col- 
lapse which is partly induced by the pressure effect. 

• As long as the external radiation is assumed to turn on when the filamentary clouds 
form, low initial density (n < 10 2 cm -3 ) filamentary clouds with moderate line mass are 
expected to fragment into very massive clumps (~ 1O 5 M ) as a result of photodissociation 
of molecular hydrogen. This result which is originally indicated in Paper I is confirmed in 
this paper using one-dimensional hydrodynamical calculations. 

• The external dissociation radiation increases fragment mass when the filamentary cloud 
fragments during adiabatic phase after the external radiation turns on. On the other 
hand, when the filamentary cloud with sufficient line mass reaches the loitering point, the 
dissociation radiation decreases fragment mass, which is consistent with Omukai & Yoshii 
(2003). 

As seen in figure 15, the thermal evolution of a filamentary cloud and a spherical cloud 
is different after each cloud reaches loitering point. The thermal evolution of the filamentary 
cloud is isothermal, and temperature of spherical cloud increases. This is explained as follows : 
since the central region of the filamentary cloud is approximately dynamical equilibrium (§3.3), 
we have 

IP GJ_ 

p r r 

T ~ fininGl oc const. (27) 

As for spherical cloud, the central region is not dynamical equilibrium. Hence, we investigate 
relation between T and n from balance between adiabatic heating and H2 cooling. Adiabatic 
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heating rate is 

-P^i^P^ocTn 1 / 2 (28) 
at p tfip 

where we assume that spherical cloud collapses in free-fall timescale. When density is larger 
than the critical density of H 2 , Ah 2 oc nT a (a ~ 3.8 at T ~ 300K and a ~ 4.8 at T ~ 1000K), 
and we have 

oc — oc V. (29) 
p n 

From equations (29) and (30), we have T oc p 1 / 2 ^ -1 ). Temperature depends on density as, 

!ni/2(a-i) sphere , , 

(30) 
const. filament. 

In this paper, we consider the filamentary cloud with variety of line mass. Hence, we 
estimate / for the filamentary clouds in cosmological simulation. As an example, we refer with 
figure 2 of Greif et al. (2008). In this figure, the filamentary cloud with density n ~ 10~ 2 cm -3 
and radius ~ 7kpc is shown. Line mass of this filamentary cloud is ~ 7.8 x 10 18 g/cm. When 
temperature is 300K, the critical line mass is / cr i t ~ 3.5 x 10 17 g/cm. Hence, / ~ 22 and according 
to our results in §4, fragment mass may increase if the external radiation turns on at n < 1cm -3 . 

For simplicity, the model and numerical calculations in this paper are one-dimensional 
for the filamentary cloud and one- zone for each fragment. In order to discuss fragmentation, 
we assume the condition for fragmentation (§2.3) and assume that each fragment is spherical. 
During further collapse, spherical clouds are possible to be the filamentary cloud again or 
may become disk-like if it rotates. Although we discussed a possible sub-fragmentation at 
the loitering point for each clumps in §4.2, the final fate of the cloud is still open question. 
Furthermore, if the filamentary cloud with large line mass is not axisymmetric, it may become 
sheet-like cloud. Such a cloud may collapse and fragment into many filamentary clouds. In 
this paper, the external radiation is assumed to be uniform, and the intensity does not depend 
on time. These problems require three-dimensional calculations. As for three-dimensional 
simulation with the external radiation, Susa (2007) investigated collapse of spherical cloud under 
the single light source. However, further investigations with three-dimensional simulation which 
statistically investigate fragment mass of the filamentary cloud under the external radiation will 
be desirable. Despite simplicity, one-dimensional hydrodynamical calculations in this paper are 
useful in the view point of extracting physical processes which are important in formation of the 
astronomical objects. These one-dimensional calculations and the realistic three-dimensional 
calculations may be complementary. 

We thank Fumio Takahara for fruitful discussion and continuous encouragement. We 
also acknowledge the referee for improving the manuscript. 
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Appendix 1. Virial equation for uniform sphere 

We multiply 47rr 3 by both hands of equation of motion, 

Dv dP GM tk , 

p m=-^-— p ' (A1) 

and integrate in respect with r. Then, the left side hand of equation (Al) is 
Dv , r Dv 



/ A-Kr 6 p-—dr = / r——dM 
J H Dt J Dt 

ID 2 



2Dt 2 

About the first term of right side hand of equation (A2) 



r 2 dM - / v 2 dM. (A2) 



Jr 2 dM = J Airr 4 pdr 



47r P o5 



5 

> 

5 



fikfi? 2 , (A3) 



1D2 V.M = I^(-M^ 



and 

2Dt 2 J ' 2Dt 2 \h' 

3 n /DR\ 2 3^D 2 R 

= -M\ +-MR -. (A4) 

5 \DtJ 5 Dt 2 K J 

About the second term, 



v 2 dM 



J Airr 2 pv 2 dr 

! il ""'\^)'h dr 



where we use the following relation, 
DR\ r 

mh: < A6 > 

since velocity is in proportion to r because of uniform density. Hence, the left side hand of 
equation (A2) is 

r Dv 3 D 2 

j^ p - ir = .„ R — R , (A7) 

On the other hand, about the right side hand of equation of motion, term of pressure gradient 
is 

f ndP , . . kaT . , . 

- / 4vrr 3 — dr = 3( 7a di - 1 — M, (A8) 
J dr /imn 
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where 7 a di is adiabatic index. The second term is 
,GM , t,. V ,G 



/ 47rr 3 — —pdr = f (47r pf : — r 4 dr 

J T J 3 



= -?^. (A9) 
5 R K J 

Finally, we have virial equation for uniform sphere, 

3 , r Dv . k B T 1 3 GM . , , 

-MR— = 3 7 adi - 1 ) — -5 ~ A1 

Dv _ 10 fc B T 1 GM 

T>i ~ Y fim R R~ ~R^' ^ ' 

where we use •y^i = 5/3. 
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